function testGc
clc
%element 1
B1 = [   -1 0 1 0 0 0;
        0 0 0 -1 0 1;
        0 -1 -1 1 1 0   ];

%test1
Gc1 = [  1 0;
        0 0;
        0 1    ];
xi1n = [1 1.2]';
d1 = [xi1n(1) xi1n(2) 0 0 0 0]';
B1*d1 + Gc1*xi1n
d1 = [-xi1n(2) xi1n(1) 0 0 0 0]';
B1*d1 + Gc1*[-xi1n(2) xi1n(1)]'

%test2
Gc1 = -[  1 0;
        0 0;
        0 1    ];
d1 = [0 0 xi1n(1) xi1n(2) xi1n(1) xi1n(2)]';
B1*d1 + Gc1*xi1n
d1 = [0 0 -xi1n(2) xi1n(1) -xi1n(2) xi1n(1)]';
B1*d1 + Gc1*[-xi1n(2) xi1n(1)]'

%element 2
B2 = [   0 0 1 0 -1 0;
        0 -1 0 0 0 1;
        -1 0 0 1 1 -1   ];

Gc2 = [  -1 0;
        0 0;
        0 -1    ];
xi2n = [1 1.2]';
d2 = [0 0 xi2n(1) xi2n(2) 0 0]';
B2*d2 + Gc2*xi2n
d2 = [0 0 -xi2n(2) xi2n(1) 0 0]';
B2*d2 + Gc2*[-xi2n(2) xi2n(1)]'

Gc2 = -[  -1 0;
        0 0;
        0 -1    ];
xi2n = [1 1.2]';
d2 = [xi2n(1) xi2n(2) 0 0 xi2n(1) xi2n(2)]';
B2*d2 + Gc2*xi2n
d2 = [-xi2n(2) xi2n(1) 0 0 -xi2n(2) xi2n(1)]';
B2*d2 + Gc2*[-xi2n(2) xi2n(1)]'
